Generation and remote detection of THz sound using semiconductor superlattices 
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The authors introduce a novel approach to study the propagation of high frequency acoustic 
phonons in which the generation and detection involves two spatially separated superlattices ~ 1 
apart. Propagating modes of frequencies up to ~ 1 THz escape from the superlattice where they 
are generated and reach the second superlattice where they are detected. The measured frequency 
spectrum reveals finite size effects, which can be accounted for by a continuum elastic model. 



With the availability of ultrafast laser pulses, there has 
been much effort in the study of the propagation of high 
frequency acoustic modes using picosecond ultrasonics [lj, 
2] . In this approach, the absorption length of the light de- 
termines the wavelength range of the acoustic waves that 
are generated. Alternatively, high frequency ultrasound 
can be excited by light pulses usingperiodic structures 
such as a superlattice (SL)@, 0, Uu 0, S, HJ- In such 
structures, the modulation of the elastic properties intro- 
duces an additional periodicity that backfolds the bulk 
phonon branches into a smaller Brillouin zone [10]. The 
new zone contains modes that have nonzero frequency 
at the zone-center and can couple to light. In contrast 
to the picosecond ultrasonics method, the wavelength of 
the waves generated is determined by the period of the 
SL. Propagating folded acoustic phonons have been ob- 
served using two-color pump-probe experiments where 
the pump and probe are absorbed in different regions in- 
side the sample Q. In this work, we demonstrate a new 
method to probe acoustic vibrations using two identi- 
cal SLs as acoustic transducers. The acoustic waves are 
generated and detected by ultrafast laser pulses in op- 
tical pump-probe experiments. We identify the modes 
observed with the aid of a continuum model, which also 
explains a fine structure that arises due to the finite size 
of the SLs. 

The time resolved studies were performed at 77 K us- 
ing a standard pump-probe setup, in which an ultrafast 
optical pulse is used to coherently excite acoustic vibra- 
tions. The probe pulses reflect off the sample, and we 
record the change in the reflectivity due to the presence 
of the coherent phonons as a function of the time delay 
between the pump and probe pulses. We use ~ 70 fs 
pulses from an optical parametric amplifier tunable in 
the visible range (400 — 650 nm). The energy of the 
pulses is ~ 80 nJ and the pulse wavelength was centered 
at 530 nm. The probe can be focused on the same side 
of the sample as the pump, or on the opposite side; see 

Fig. rjj 

We studied two samples grown by molecular beam epi- 
taxy on a (100) GaAs substrate that consist of a thick 



GaAs layer of 1.2 (im or 0.6 /im between two identical 
superlattices of 25 periods of 12 A-thick GaAs/34 A- 
thick AlAs. In what follows, we refer to these samples as 
sample A and B, respectively. To access the second SL 
we removed the substrate by mechanical polishing and 
chemical etching. Our scheme makes use of the photoc- 
lastic coupling of acoustic vibrations to light through the 
zone folding of the acoustic branches. The absorption 
of the pump pulse in one of the SLs launches a high fre- 
quency acoustic wave that propagates towards the second 
SL where it is detected by the delayed probe pulse (Fig. 
[lj. This particular design makes it possible to generate 
and detect vibrations in spatially separated regions. 

The samples were characterized by Raman spec- 
troscopy. Spectra obtained with the 514.5 nm line of 
an Argon-ion laser in the backscattering geometry show 
a doublet near ~ 1 THz which corresponds to the first 
zone-center folded acoustic modes [ll|. The frequencies 
obtained from these experiments match the calculated 
ones using an infinite continuum elastic model, originally 
developed by Rytov [12j . 

In Fig. [2] we show experimental differential reflectivity 
data for the two samples. In (a) we show a typical time 
trace for the configuration in which the pump and the 
probe arrive at the same face of the sample. The data 
clearly exhibit a high frequency oscillation superimposed 
on top of a large, lower frequency vibration due to stimu- 
lated Brillouin scattering [1(. Results for the situation in 
which the pump and the probe strike on opposite sides 
of our structure are shown in Fig. [2] (b). For samples A 
and B we observe, respectively, an oscillating transient 
in the reflectivity arriving after time delays of ~ 115 ps 
and ~ 230 ps. These values agree well with the time-of- 
flight of acoustic waves traveling across the bulk layer. 
We note that the high frequency folded phonons arrive 
at the second SL after the transit time, as shown in the 
inset of Fig. [2] (b), which is an expanded view of the 
region indicated by the dotted box. The amplitude of 
the oscillation due to the folded phonons is nearly the 
same in the data of Fig. [2] (a) and Fig. [2] (b), indicating 
that they propagate without a measurable attenuation 
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through the GaAs layer. 

Figure [3J (a) shows the Fourier transform of the trace 
presented in the inset of Fig. [2] (b) . There are four main 
features present .a large amplitude vibration due to Bril- 
louin scattering [l( at 70 GHz (not shown in the figure) 
and a triplet due to zone-center folded phonons with com- 
ponents near 1 THz. Two modes of the triplet correspond 
to the back-scattering wavevector (BS) while the remain- 
ing one (FS) originates from the zone center mode that 
is Raman active The frequency of the BS modes is 
determined by the wavelength of the probe pulse. This 
dependence can be understood from the following argu- 
ment. The change in reflectivity can be written as [l[ 

/+oo 
e 2ik " z P(z)r](z,t)dz, (1) 
-oo 

where rj(z, t) — du(z, t)/dz is the time dependent strain, 
P(z) is the corresponding component of the photoelas- 
tic tensor and k p is the wavevector of the probe. In 
an infinite structure, Eq. ([T]) has non-zero contributions 
from those components of r\ that satisfy the phase match- 
ing condition q = ±2k p . This means that, regardless of 
the generation mechanism, the probe selects only those 
modes that satisfy wavevector conservation. Further- 
more, the propagating strain pulse moves at the speed 
of sound and, therefore, the corresponding frequency for 
the phase matched low frequency oscillations in Fig. [2] 
is z^br = kp/nv where v is the speed of sound of the 
SL. For the same reasons, the folded phonon doublet at 
~ 1 THz is composed of modes which, in the extended 
zone scheme, have wavevectors q — 2t:/D ± 2k p , where 
D is the SL period. 

We calculated the change in reflectivity induced by the 
pump-generated acoustic vibrations. The corresponding 
equation of motion is [l[ 

d 2 u(z,t) da(z,t) 

pM-ap- = -aF- (2) 

where p is the density, <j{z,t) = C(z)du(z,t)/dz is the 
stress, C(z) is the elastic stiffness, and u(z,t) is the 
atomic displacement. The energy deposited by the laser 
pulse sets up an initial stress (To (2) = K(z)\E(z)\ 2 , where 
E(z) is the electric field of the pump pulse and K(z) is 
a square wave that represents the absorption in different 
layers of the SL *2]. The use of a continuum model for 
structures with such small periods is somewhat arguable 
since the bulk properties of each constituent are not well 
defined and, hence, the local index of refraction or, equiv- 
alently, the absorbed energy in each layer is not known. 
With this in mind, we adjusted the ratio ATcaAs/^AiAs to 
reproduce the relative intensities of the observed modes. 
Regardless of the generation mechanism, the waves are 
detected by stimulated light scattering through the pho- 
toelastic mechanism, as stated in Eq. (J]). 

Calculated spectra are shown in Figs. [3] (b) and (c) 
together with the experimental result (a). The best fit 



was obtained for i^GaAs/^AiAs = 1-2. For clarity, we 
have included the dispersion relation for the infinite SL 
(lower panel) . In the spectrum corresponding to a single 
semi-infinite SL, Fig. [3] (c), the flat features with zero 
intensity at 0.5 THz and 1 THz reflect the acoustic mini- 
gaps whereas the doublet at ~ 1 THz corresponds to 
the folded phonons observed in back-scattering Raman 
spectra (lp| . 

Figure [3J (b) shows the calculated spectrum for the fi- 
nite SL in the case where the probe is incident on the 
second SL. The calculations exhibit a complex structure 
with an overall modulation of the intensity, which is also 
visible on the lower side of the measured spectrum in 
(a) . The frequency spacing of the comb-like structure is 
Av ~ 20 GHz. This value corresponds to the inverse 
of the transit time of the acoustic pulse through the SL. 
This periodic modulation manifests itself in the time do- 
main traces as a time delayed signal. We note that the 
high frequency oscillations shown in the inset of Fig. [2] 
are delayed from the main Brillouin oscillations. This de- 
lay, of t ~ 50 ps, is consistent with the frequency spacing 
of the comb. 

The folded phonons in Fig. [3J (b) also show a fine 
structure which arises due to the finite size of the SL. In 
the experimental data, this fine structure is particularly 
evident in the higher frequency BS-peak. In addition, 
the peak at the zone center (FS) is also visible. This 
feature appears stronger in the experiment than in the 
calculations. We tentatively attribute this discrepancy 
to the interface contribution to Ar[9(, which is not taken 
into account in Eq. (TT]). 

In conclusion, we introduced a new method to probe 
the dynamics of acoustic phonons of frequencies in the 
THz range using two spatially separated SLs. The fre- 
quency spectrum of the detected waves has a complex 
structure that originates in the finite number of periods 
of the SLs. A continuum elastic model accurately pre- 
dicts the observed frequencies as well as the fine struc- 
ture. However, some parameters in the model had to be 
adjusted to obtain the correct amplitudes in the spec- 
trum. 
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FIG. 1: Schematic diagram of the double SL structure and 
the geometry of the pump-probe experiments. The samples 
consist of a thick GaAs layer of 1.2 jj,m or 0.6 /xm between 
two identical SLs. 



FIG. 2: Time-domain data for samples A and B. Pump 
pulses generate sound waves on the front superlattice, that 
are later detected on the front (a) and the backside (b) of 
the sample. Trace (a) and the inset of (b), obtained on sam- 
ple A, show folded-phonon oscillations superimposed on the 
larger, low-frequency oscillations associated with stimulated 
Brillouin scattering. 
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FIG. 3: (a) Fourier transform of the time trace for sample A. 
We observe three peaks that correspond to the backscattering 
(BS) and forward scattering (FS) modes accessible in Raman 
experiments; (b) calculated spectrum for the double SL struc- 
ture, where the probe is incident on the second SL and (c) 
calculated spectrum for a semi-infinite structure showing the 
BS doublet. The lower panel shows the dispersion relation for 
an infinite SL. 




Front SL GaAs 




Back SL 




delay [ps] 




delay [ps] 




Frequency [THz] 



